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Abstract 

We consider the critical alternating quantum spin chain with ^ spins. Using 
the Bethe ansatz technique we find explicit expressions for the 5'-matrix of the model. 
We show that in the limit that g-t — > oo our results coincide with the ones obtained 
for the principal chiral model level one, for the LL (RR) LR scattering. We also 
study the scattering of the bound states of the model and we recover the results of 
the XXZ (sine -Gordon) model. 
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1 Introduction 



In the famous paper [0] the principal chiral model (PCM) with SU{2) symmetry with 
Lagrangian density 

Lo = ^tr{d,g-\x)d''g{x)) (1) 
were firstly shown to be equivalent to a Gross-Neveu (GN) fermionic model 

N _ ( ^ - \ ^ 

a=l \a=l ) 

in the limit of infinite number of fermionic colors (r^are Pauli matrices here). Secondly, 
in order to avoid ambiguities related to the axial anomaly, the fermionic model was con- 
jectured to be equivalent to the single fermionic field with isotopic spin S {N = 2S). The 
lattice analogue of the last model is the spin S quantum spin chain exactly solvable by the 
Bethe ansatz method. This scheme was successfully applied to PCM and GN models with 
different symmetries 0, 0, |Q, 0. Another modification of this scheme can be used for 
analysis of the PCM with topological term 

Li = ^tr{d^g-\x)d^g{x)) + ne,, f dtdPx tr ig-^gg-'d^ gg'^d" g) (3) 
zckq Jo ^ ^ 

where in the second term g{x^ t) interpolates between g{x^ 0) = 1 and g{x, 1) = g{x) and n 
is an integer. This model was shown to be equivalent to fermionic GN model in the limit 
of infinite number of fermions, with broken chiral symmetry, i.e. different number of colors 
for left and right fermions: 

L, = f: uid^ua + E vld-va + Ao E E ("I^^^^) (^^V,) (4) 

a=l a=l a=l 6=1 

where Ua{va) are the left (right) Dirac fermions, Aq = A* - the generators of symmetry 

group and the difference of fermionic color a numbers is q+ — q- = n. Here we are focused 
on the case where g+ — g_ = 1, which corresponds to the PCM level one. The corresponding 
statistical model for this case is the critical alternating spin chain with ^ spins (see 
e.g. [0). We are going to find the explicit S'-matrix for the spin chain, using the Bethe 
ansatz method, and we will reproduce the proposed massless S'-matrix for the model (0), 
see i. 
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2 Bethe ansatz equations 



2.1 Ground state and excitations 

We focus here on the alternating spin chain with ^ spins. In order to introduce a 
mass scale to our system we consider the chain with inhomogeneities [@, 0, namely 

^. = (-r-, (5) 

where j is q± and let us also consider that g+ is even (obviously g_ is odd). In other 
words in the sites of the chain with ^ spin Ug^ = ^ whereas in the sites with ^ spin 



Uq_ = — Following the standard Bethe ansatz technique 0, [§-[|14| for the model (0), 
one gets Bethe equations in the form 

1 1 *^ 

e,4A, - -)^+e,.(A« + = - H e^iXa - Xp) (6) 

where 

sinhMA + f) 
^"^^''^^^ sinhMA-f) 

where = ^ is the anisotropy parameter. The spin, energy and momentum of a state are 



characterized by the set of quasi particles with rapidities (BA roots) Aq, 0, |]T2|, [|T5 



rn sin /in 



2vr ,=1 sinh/i(A, - ^ + f ) smhfi{Xj - LI 

M o^riV, n(\. - Dl 1 ^" 



2 > 



p = -y y log — ^ — ^^2^ — ^ (9) 

2^Un4:,,- %inhMA,-^-f) 

5^ = N^q_ + Ar_^g+ - M. (10) 

We consider here the special case where = = A^. We are interested in the limit 
N±iQ±ii^ ~^ oo. The thermodynamic limit N± — > oo of equation (^) can be studied by 
standard methods. The string hypothesis says that solutions of (^ in the thermodynamic 
limit are grouped into strings of length n with the same real part and equidistant imaginary 
parts 

X(^^j) = Xl+'-{n + l-2j) + z^{l-v), j = l,2,...,n (11) 
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where is real and v — ±1 determines the parity of the string (positive negative parity 
string). We assume that v > q±. The ground state consists of two Dirac "filled seas" by 
strings of length q±. Therefore the Bethe ansatz equations for the ground state become, 

n x^AK-^-^f= n n^n.(As-A^^) (12) 

3=1+, Q- j=q+,<i- /3=i 

where n can be q±, and 

n—m+l\ 

(A)e 

|n— m+3| 

(A)... (A) 6(n+m— 1) (A) 

-E^nm(A) = e|.„_„|(A)ej^„_„^2l(A) • • • e(„+m-2)(A)e(n+m)(A). (13) 

For what follows it is necessary to introduce the following notations 

/\ N /X I cosh//(A+ I) 

3„(A;.)=e„(A±-) = ^^^j^^^, (14) 

GnmW = g\n-m\ W9'\n-m+2\ (A) ■ • • fi'^n+m-2) (A)^(n+m) (A) , (15) 

a„(A;i/) = ^^iloge„(A;i/), 6„(A; z^) = ^^ilog^„(A; i/) (16) 

and the Fourier transforms of a„, 6„ are given by 

smh((u - n)^] 

«n(c^; ^) = < n < 2i., (17) 

smh(-|^j 



sinh(^f ) _ sinh ((n - 2«y 



'"^"'"^ = -ihn#) 0<-<-'= sinh(f) -<-<3- (18) 

We also need the following expressions 

(Z„^(A), Anm{X), 5„^(A)) = ^^ilog(X„^(A), Km(A), G'nm(A)), (19) 

where the Fourier transforms of the last expressions are 

sinhf (i/ — max(n, m))%] sinhf (min(n, m))%] 

= — lh(f)ih(i) ■ 

2 coth(^) sinh(^(z^ — max(n, m))^^ sinh(^(min(n, 77i))|) 
^nm('^) = sinh(^) ^^■^^ 



3 



2 coth(^) smh(^) sinhf 

^""■<-) ■ 

Finally, the energy (|^) takes the form 

E = - E Y.Z^AK-^-^)■ (23) 

Now we are going to solve the Bethe ansatz equations. We take the logarithm and the 
derivative of the Bethe ansatz equations and we also use the Maclaurin expansion 

E /(A,) ~ N / /(A)a(A)dA - ^ /(A,), (24) 

where cr is the density of the state and u is the number of holes. Then for the densities of 
the ground state (no holes) we obtain the following integral equations, 



<W= E ^n,(A-^)- E (^n,*a^)(A) (25) 
where * stands for the convolution and n can be Qj^ or q^. We can easily solve the last 



a-{\) = six - i-^) = ^ , , , (26) 



equations (25) and we find that 

ao 2coshf7r(A 
where its Fourier transform is 

s{u)= E (4A,)(u;)= (27) 
Here i? is the inverse of the kernel K of the system of the linear equations (pSf ), in particular, 

^nm(^) = (1 + Am(t^))(^nm + ^nm(w)(l " f^nm) (28) 
Rnm{uj) = E ((1 + -4jj(a;))5„m(l - (^nj) - ^nm(t^)(l - 5nm)), (29) 

where the determinant of K is, in terms of trigonometric functions, 

4coth2(^) sinhf(z/ - q+)^) sinh(g_^) sinh(^) 

— ^-^mW) ' -■ 

Now we consider the low lying excitations of the model. There can be holes in both seas of 



q± strings and also strings of other length and parity (see equation (pJ])). Let us consider 



the state with u± holes in the q± sea. The contribution of the holes modifies the densities 
of the states in the following way (see (plD), 



a"(A) = e"(A) + 1 ^ E irr(A - A^), (31) 



where e"'(A) = s(A — 



(-)" 

ao 



N . 

i-=<i+,<i-' "=i 



Kr{^)= E (32) 



In particular, K in ( pT]) have the following explicit form in terms of trigonometric functions 



sinh (i>-2)^ . 1 

' ^ ^ 2cosh(|)sinh((/>-l)|) ' ^ ' 2cosh(f) ^ ^ 

with z/ = i> — q_, and also 

1 ^ sinhf(g+ -2)^) 

^ ^ 2cosh(f)' ' ^ ^ 2cosh(f)sinh((g+-l)f) ^ ' 

The energy of the state with u± holes in the q± seas is given (see (p3D) by 



^ = ^0+ E Ee'(A'j, (35) 
j=q+,g- a=i 

where E'o is the energy of the ground state and e^(A) is the energy of the hole in q± sea 
(we write e"*^ instead of e''± for simplicity). Finally, we compute the spin of the holes, and 
we can see that the spin of a hole in g+ sea is from (10) s"*" = ^, where / is just a factor 



/ = 1 + -S^ which is obviously 1 for z> oo. Therefore, we consider the spin of the hole 
to be I for what follows. The spin of a hole in sea is see (|10D s~ = 0. Finally, we 
conclude that the hole in the q± sea is a particle like excitation with energy e^, momentum 

= n (-)"^^ ' P^^^^ = -J + ^tan-^(sinh7r(A - ^)) , (36) 

2cosh(^7r(A - i^)j 4 2 ^ ao ^ 



and spin = |, s =0. We can easily check that in the scaling limit where A << ^ the 
energy and momentum become from (|^) and (|36[) , see also |@|, 0, 



e+(A) =p+(A) ~ e"^e"\ e-(A) = -p-(A) ~ e"^e-^\ (37) 



These are the energy and momentum of the "right" and "left" movers respectively (see 
e.g. 0). The factor e °o provides a mass scale for the system. 

As we mentioned before there are also states that consist of strings. There exist two 
types of strings as we see from equation ( [TT| ) (positive and negative parity). The presence 
of the string gives rise to some extra terms in the Bethe equations which give the 1/N 
contribution to the density, in particular the Bethe ansatz equations in the presence of a 
positive parity /-string become, 



n x„,(A^ - i-^)^ = n ni?„,(A^-A^,)E„KAs-Ao) m 

j=g+,<i- ° j=g+,g- /3=i 

where Aq is the real center of the string. The densities of the state are 

a"(A)= E ^n,(A-^)- E (^n.*^')(A)-^^nKA-Ao), (39) 

Cin • . iV 

J=Q+,Q- " 3=-,+ 

we solve the system of the linear equations and we find that: 

a"(A) = e"(A)-lK^'(A-Ao), (40) 

where 

kl'{u)= Y: (4-^n,)(^). (41) 

j=Q+,g- 

We have to consider two different cases for the positive parity strings. In particular, we 
choose the following strings, namely, the / = g+ + 1 and / = g_ — 1. For each string we 
obtain different densities, i.e. for the / = g+ + 1 string, 

kl-^-'-^'\u) = = (42) 

smh(^(z^ - 1)1 1 



and for / = (/_ — 1 

. t u f ^^ sinhf(g+-2)^ 

sinh((g+-l)f)) 

The energy of the state from (^, (^Tj) and (P^), is given by, 

E = - E E%(AL-^)- E %(Ao-^), (44) 



we find that E = Eq, therefore we conclude that this is a state with zero energy contribu- 
tion. The spin ([ToD of the + 1 string is s = — 1 (again we omit the factor / for the spin) 
and the spin of the q- — 1 string is zero. Finally, we study the state with an one negative 
parity string. The Bethe ansatz equations then become 

n x„,(A^-^r= n n^n.(^2-A^^)G„i(A^-Ao) (45) 

Similarly to the previous case the densities of the state are 

a"(A) = e"(A)-i^ir!!i(A-Ao), (46) 

where 

k^\u)= Yl (A,^n,)M- (47) 

j=g+,<i- 

Then the contribution of the densities becomes 

kl^\u;)= '^^^^^ kl-\iu)=0. (4J 

sinhf (z/ — 11 '" 
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The energy for the string is given by 

E = - E E%(^«-^)-^.+ (Ao--)-V(Ao + -), (49) 

i,j=q+,q.a=l ^0 



and one can show from (|49|) , (^Gf ) and (^) that the negative string has also zero energy 
contribution. Again the spin of this string is s = —1. These are all the necessary states 
to describe the scattering processes for the specific model. Note that the inhomogeneities 
modify only the energy and the momentum of the model. However, the contributions 
to the densities, which essentially determine the scattering amplitudes, remain unaffected 
by the presence of the inhomogeneities. In other words the computation of the S'-matrix 
is not affected by the presence or not of inhomogeneities in the spin chain. Therefore, the 
iS-matrix holds true for any A, and we do not have to consider the scaling limit. 

2.2 Physical S matrix 

Having studied the excitations of the model in the previous section we are ready to compute 
the complete iS-matrix. To do so we follow the formulation developed by Korepin, and later 
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by Andrei and Destri [16|. First we have to consider the so called quantization condition. 

(e'^^PS -l)\XuX2) = (50) 

where p is the momentum of the particle, the hole in our case. For the case of two holes 
in g+ sea (this is the triplet state s = 1), we compare the integrated density (|31|) with the 
quantization condition. Having also in mind that, 

e"(A) = p-{X) (51) 
we end up with the following expression for the scattering amplitude 

5o(A) =exp{27r2^ / Kl+^+ {\ - \.i)d\] (52) 
i=i 

where A = Ai — A2 and Ai,A2 are the rapidities of the holes. In terms of trigonometric 
functions (|33D 



/ r -.^A sinh((z>-2)|) dcu 

So{X)=exp\- e ^ 7 ^ — 

^ J-00 2coshf^)sinh(rz>- 1)^) ^ 



2cosh(f)sinh(^(z>- l)f' 

00 r('JA_ I 2.7-2 I 1^p('_^A_ I _2j_\ I 2j-l I 1 2j-l N 

,=1 r(4\ + ^ + i)r(4\ + ^) r(^ + ^ + i)r(^ + ^)' 

To get the singlet we consider the state with two holes in the g+ sea and the + 1 string. 
Looking at the spins of the hole in g+ sea and the spin of the g+ + 1 string we conclude 
that this is indeed a state with spin zero. Note that the real center of the string is in the 
middle of the two holes (Aq = ^^). We conclude from (^, (gS]) and (|0|), (g^) that 

5+(A) = exp{-2^2 K^+^^++')(A - Ao)rfA}5o(A) (54) 
and after some algebra we get, 

sinh(27fTT(A + «)) . 
sinh(^2{s3i)(A-^)j 

However, we need another state of spin zero see [0, and this is the state with two g+ holes 
and one negative parity string with its real center in the middle of the holes. For which 
we obtain (EGD-fHHD 



coshf (A + 
S.{X) = X.r^ A SoiX). (56) 
cosh(2(^(A- 



The equations (|53[), (|55D and (0) give the XXZ S'-matrix (sine -Gordon S'-matrix, provided 
that 0^ = 8n(l — 4) and A = ^) ||], ]12|, |18[. Moreover, we have holes in the g_ sea, 



which are particles with zero spin. The scattering between two such holes (see (0), (M)) 
is described by the following S'-matrix 

r r°° ■ X sinhf (g. — 2)^ ) duj. 

S',{\) = exp|- / e-^-' , V — } , (57) 

^ -^-oo 2cosh(f)sinh((g+-l)f) ^ J 

which, in the limit q± — > oo, corresponds to the XXX scattering amplitude []^, namely, 

S'CX) - expl- r e--^:^^) - (58) 
6o(A)-exp| 3 cosh(f ) u; / " r(:#)r(i + f ) ' ^''^ 

If we add a g+ + 1 string in the above state the only thing that changes is the spin of the 
state which becomes -1, the scattering amplitude remains unaffected. There is one more 
state which completes the correspondence to the XXX S'-matrix: the state with two holes 
in the g_ sea and a positive parity g_ — 1 string (spin zero state). For such state the change 
in cr'^- density gives rise to an extra phase, namely, 

S'4X) = ^S',{\). (59) 

The last two equations give the XXX S'-matrix (see e.g 0, 0). The equations (|53D - (|59D 
give actually the correct Sxxz x Sxxx matrix [^. In the limit v ^ 00 the XXZ matrix 
reduces to the XXX S'-matrix, i.e. the matrix given by (|55|), (p^. Therefore we end up 
with the SU{2)l x SU{2)r S-matrix [0. Finally, let us consider the state that 

consists of two holes in g+ sea and two holes in g_ sea (we consider even number of holes 
in each sea). The scattering between the two different holes is described (^I]), (^3]), 



~ r ~ 1 dill ^ 

5Lii(A)=exp{-/ e-^— —-- — }, (60) 
j-00 2cosh(|) oj > 

which can be written in the following form 

Sz.«(A)=tanh^(A-^), (61) 

the same result (up to an overall factor) is obtained for the spin 1/2, 1 alternating isotropic 
spin chain see [|l^. The last result (|6TD agrees with the proposed S'-matrix for the left right 
scattering [^. We identified the SU{2)l x SU{2)r S'-matrix and the Slr matrix as well. 
This is the proposed S'-matrix that describes the massless scattering for the PCM level 
one. 
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3 Bound states 



It is well known for the XXZ model that for specific values of the anisotropy parameter 
there exist bound states of "kink-anti kink". The same happens for the alternating spin 
chain as well. In particular, for the Sxxz matrix (|53|), (^5|), (|56|) and for g+ + l > ly, {u < 2 
which corresponds to < Att) there exist poles in the physical strip, Xp = i — ik{i) — 1), 
A; = 1,2, ... < -J^^ see also e.g. [1^ for a more detailed analysis. We are going to study 
here the first breather k = 1 which corresponds to a negative parity string of length one. 
For i> < 2 (obviously here we keep i) finite) the Fourier transform of Biq^ changes, namely 

2 sinhffzy - coshffi/ - 1)^) 

- — ^ — siffr — 

Having that in mind and also the equations (^61), (0) we conclude that the jj contribution 
of the densities for the breather state are, 



. , cosh (/>-l)f . , cosh (^-2)^ 

^ ^ cosh(f) ^ ^ cosh(f) ^ ^ 

The negative parity one string for z/ < 2 is a state with positive energy and zero spin and 
it corresponds to a breather. In particular the energy of a breather is see (^), (0) 

= ^(co.h(.(A - + ^2))) co*(.(A - ifca- ^ • 

where 



cosh (i>- 2)^, 
cosh(l^) 

it is indeed positive and it is the same as in the case of the critical XXZ model in the 
"attractive regime" (see e.g. |12[, pO[). In the scaling limit A << ^ the energy and the 
momentum of the first breather become, 

ei(A) =Pi(A) ~2sin(^(z>-l))e-^e-\ (66) 

this is a dispersion relation for massless particle where the factor 2 sin(|(i> — l))e °o is the 
corresponding mass scale. 

^In the XXZ critical model with inhomogeneities uij — (— )-'6 (see e.g.|9|, [||) one can show that the 
mass of the first breather is exactly 2sin(^(i> — l))2e^'^®, which is also the mass of the sine-Gordon first 
breather, provided that 2e^^® is the "kink" mass. 
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We can write the Bethe ansatz equations for a state with two breathers, 

n 9A^--^f= n nG'i,(A-Ai)ne2(A-A.). (67) 

j=g+,g- ^ j=g+,Q- a=i i=i 

The density will be 

-a{X) = b,^{X-^) + b,_i\ + ^)- Y: (5i,*aO(A)-i-^a2(A-A.) (68) 

«0 "0 j=q+,q- i=l 

which we can write in the form, 



&iX) = e,iX) + ^j:K,iX-X,). (69) 

1=1 

The fourier transform of Big is given by equation (p2D, whereas the Fourier transform of 
Biq^ is given by (0). Therefore the Fourier transform of the Kf,, after some algebra, is 



cosh (2^-3)^ 

^^'(^> - — 

To study the scattering between two breathers we consider the following quantization 
condition for a breather 

(e2^^^^5-l)|Ai,A2) = (71) 

where pi is the momentum of the breather, then if we compare the last equation with 
(^) (after we integrate it) we conclude that 

5.(A) = exp{ r e--f!!K(^£zM!!^} = sinhM)+.sm(.(g-l)) 

U-oo cosh(|) u; J smh(7rA) - zsm(7r(i/ - 1)) ^ ' 

The result agrees with the one obtained in XXZ (sine-Gordon) model, see e.g. |12|, |1T8| , 



]. Here, we restricted ourselves to the case of the "elementary" breather only. The 
above analysis can be generalized, in a similar way, for any general n-breather. 

4 Discussion 

We considered the alternating critical quantum spin chain with ^ spins. We studied 
the excitations of the model and their scattering. We showed that in the limit that q± ^ oo 
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the 5'-matrix we found coincides with the proposed one for the PCM level one. The 
natural next step is the generalization of the above method for PCM model at any level. 
In this case, one expects a relevance of alternating fused RSOS models by analogy to |jl9 . 



A Bethe ansatz analysis similar to the one presented here, was done by Wiegmann for 
the reproduction of massive matrix for 0(3) a model with ^ = 0. It was done on the 
basis of special parameters limit of the spin S XXZ chain. In this context the interesting 
problem is a relation of the alternating spin chain to the massless S'-matrix of 0(3) a 
model with ^ = tt |^ and to its relatives, recently proposed in for other symmetric 
space (T-models with 9 term. It would be also interesting to study the thermodynamics 
of generic alternating spin chain (the case of spin 1/2 spin 1 chain was studied in p3| . 



2J], p5[) in order to see its relation to the thermodynamics of the corresponding field 
theoretical model, which was studied by methods of TBA. We hope to report on these 
issues in a future work. 
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